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O , Abstract. E x B-drifting jets have been generally ignored for the past 25 years even though they 

'. may well describe all the astrophysical jet sources, both on galactic and stellar scales. Here we present 

a \ closed-form solutions for their joint field-and-particle distribution, argue that the observed jets are near 

equipartition, with extremely relativistic, mono-energetic e^-pairs of bulk Lorentz factor 7 ^ 10'', and 

^ ■ are first-order stable. We describe plausible mechanisms for the jets' (i) formation, (ii) propagation, 

and (iii) termination. Wherever a beam meets with resistance, its frozen-in Poynting flux transforms the 

00 I via single-step falls through the huge convected potential. 



delta-shaped energy distribution of the pairs into an almost white power law, E Ne ~ E ^ with e ^ 0, 



m 
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(3 ! 1. Background 

^: 

Ph| Pair-plasma jets with ultra-relativistic bulk motion have been proposed twenty-five years ago by one 

P ■ of us (Kundt, 1979), then jointly elaborated by us (Kundt & Gopal-Krishna, 1980). They were also 
proposed by Morrison (1981), but have usually not been mentioned (cf. Begelman et al, 1984, 1994). 

c^ . Despite intermediate progress reported, e.g., in Kundt & Gopal-Krishna (1986), Blome & Kundt (1988), 
Baumann (1993), and in Kundt (1996, 2004), they have been treated with a healthy scepticism by the 
k3( \ community, cf. Blandford (2001), and Beresnyak et al (2003), with a few notable exceptions, e.g. 

5^ ■ Reipurth & Heathcote (1993), Scheuer (1996), Prieto et al (2002), Brunetti (2002), and Stawarz (2003). 
A possible reason for this lack of widespread acceptance may have been a concern about whether or not 
the beams allowed a stable transport of a broad energy distribution of high-energy charges, in the form 
of an ordered E x B-drift. The beams are indeed unlikely to transport a broad distribution. 

Instead, their Poynting-flux-flooded formation regions are expected to generate particle distributions 
at least as sharp in 4-momentum as relativistic Maxwellians, and an onsetting E x B-drift will further 
sharpen the narrow distribution towards a delta-type one. Such equipartition pair-plasma flows convect 
half of their energy as a stationary Poynting flux which is ready - wherever stalled - to broaden the 
particle distribution into an almost white power law, starting with Lorentz factors of order 10^ at their 
bottom end, and extending up to Lorentz factors of order 10^, in the form of a long high-energy tail 



whose radiated power peaks at the top end whereas its energy density peaks near its bottom end. During 
undisturbed propagation, such beams are loss-free on Mpc scales except for minor inverse-Compton 
losses on the radiation background. Note that peripheral tapping of a beam can reveal its monoenergetic 
distribution, as in Sgr A*: the convected fields vanish at the channel wall. 

We shall present exact solutions for such monoenergetic beams in section 2, and show that they 
are stable to first order. In the two subsequent sections, we shall offer reasonings why such beams are 
expected to form naturally around magnetized rotators, i.e. around rapidly rotating stars as well as in 
the centers of galactic disks, and why their radiation is expected to take the form of broad power-laws, 
from the radio to the X-ray and gamma-ray regime, with certain emission dips and excesses which can 
be understood as due to strong anisotropics in their emission patterns at high frequencies. This uniform 
model covers the observed jets from (a) newly forming stars (or YSOs), (b) forming white dwarfs, inside 
planetary nebulae (PNe), (c) young binary neutron stars (within light or heavy accretion disks), and from 
(d) the nuclear-burning centers of galactic disks (or AGN). 

Note that monoenergetic relativistic electron beams have just been produced in the lab: Katsouleas 
(2004). A few extreme and/or controversial jet sources are discussed in section 5. 

2. Solving the Beam Equations 

As has already been argued - and will be elaborated in the next two sections - the beams of the jet 
sources are expected to consist of overall electrically neutral and current-free configurations of electrons 
and positrons at large Lorentz factors 7 > 10^, convecting toroidal magnetic fields and 'radial' electric 
Hall fields w.r.t. the roughly cylindrical geometry of a beam segment. In reality, such segments of 
stationary flow have an approximately conical shape, but will be approximated by us, for simplicity of 
presentation, by cylinder segments. Note that even a conical beam does not have (adiabatic expansion) 
losses when moving through a (strictly) vacuum channel. 

We describe a beam segment by cylindrical coordinates z, s, (/?, with z growing parallel to the beam 
axis, and s, ip being polar coordinates in the cross-sectional planes. In these coordinates. Maxwell's 
(stationary, axially symmetric) equations V ■ E = Aixp , V x B = (47r/c)j ,VxE = = V-B yield 
respectively: 

ds{sE,) = 4nsp (1) 

ds{sB^) = 47rsp(3, , dsB, = -Anpf]^ (2) 

with /3 := v/c , and with all quantities only depending on the radial coordinate s. We restrict our attention 
to E X B-drifting charges, for which the acceleration 

c(7/3)- = (e/m,)(E + /3xB) (3) 



vanishes, yielding 

Es = B^(3, - B,P^ . (4) 

For realistic solutions with net charge and current zero, we have the additional boundary conditions 
at s = and R (:= beam radius): 

{sE^m = = {sEs){R) , 

(5) 
isB^)iO) = = isB^)iR) , 

and the general stationary, cylindrically symmetric solution can be Fourier-expanded w.r.t. s/R: 

sEs, sB^ ~ J]k>iCksin{k7rs/R) , (6) 

of which the first (ground) term already contains most of the physical information. 

For large 7 (>10^), we have |/5^| << /^^ ^ 1, and the leading term of the expansion reads: 

Eg ^ B^ ^ C sin{TT s / R) / s , p ^ {irC / R)cos{tts / R) / s ^ j-Jc , B-^ ~ const . (7) 

It expresses a uniform flow in z-direction, with a positive {or negative} net charge density inside of 
R/2, and the opposite charge density dominating for s > R/2 , and correspondingly with a net {positive, 
negative} current density {inside, outside} of s = R/2 , both of which peak at both ends, on the axis as 
well as at the periphery; see Fig. 1. 

This solution is as simple and transparent as one could have imagined; the involved charge densities p 
are a tiny fraction of those composing the beam: p/erif. = 10~^'^73/L44 , i.e. correspond to a minute dis- 
tortion of charge neutrality. For a jet of power L and cross-sectional area A , (ram) pressure equipartition 
(among particles and fields) requires 

7nemec2 = {E^ + B^)/8n ^ B^/An ^ L/Ac (8) 

where Ug is the electron number density. (Equipartition is plausible from what we have said, and conforms 
with the observations, cf. Begelman et al, 1984). Note that the convected electric potential $ = J Egds 
for typical jets can be gigantic, of order e$ ^ e(7rL/c)^/^ = 10^^'^eV L44 ; we shall see below that it 
can generate a high-energy tail reaching up to electron Lorentz factors of 10^ and more. Note also that $ 
is scale-invariant: A has dropped out, via s^; $ depends solely on the source power L, which is still large 
for stellar jet sources, with L ^ lO^^^^erg/s. 

Realistic jets have /S^ < 1 and P^ ^ 0. For them, a straight-forward calculation, starting from equ. 
(4) and using the integrated equs. (1,2), leads to 
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Figure 1 : Cross section through (the ground mode of) a beam segment, showing the radial dependences 



P/?^ = d,{l dss-''d,{\ dspsf - dsiJdsp/S^sfY/^ ^ dslJdss-^dMdspsf]^/^^ , 



(9) 



the latter for /J^ ^ 1 — 1/27^. It shows that /3<^ is small of order I/7 , i.e. that ordered spiralling of the 
charges should be unimportant in high-energy jets. This result must not be confused with the existence 
of helical beams, in interaction with turbulent (heavy) environs. 

So far, we have assumed strictly monoenergetic beams, with (3^ = const, which cannot be expected 
under realistic conditions. Charges whose Pz deviates from (the local value of) Es/ B^ will violate equ. 
(4), and start moving radially inward or outward, depending on their sign, whereby both {hard, soft} 
charges move {inward, outward} for opposite signs. In each case, a glance at equ. (3) shows that the radi- 
ally moving charges of deviant 7 must fall through the electric potential $ such that their 4-momentum is 
adjusted to the (locally) appropriate value for a stationary drift: more energetic ones lose, less energetic 
ones gain in energy, independent of the sign of their charge. This key stability is intuitive already on 
energetic grounds, from the shape of $, but follows directly from (3) because any radial drift implies an 
acceleration in iz-direction, via the e(E + /3 x B)-term (whose sign changes with the sign of e, and 
likewise with the sign of /3 — < /3 >). 

Note that in principle, jSz could have been radius-dependent so that a finite spread in 7-s is transported 
by the beam. But such a fine-tuning of fields and particles is unlikely to be stable, after what we have just 
found: The wings of a distribution (in 7) are removed during short distances via radial falls (of proper 
sign) through the convected E-fields - short of order 10~^° R because the E-fields are gigantic, as shown 
below equ. (8) - so that individual particle energies are stabilized quickly. 



3. Forming the Jets 

We shall now argue that monoenergetic beams - as have just been considered - are expected for the 
cosmic jet sources, i.e. are not grossly overidealized. 

To begin with, there has to be an abundant source of relativistic e^ : Magnetic reconnections, leading 
to cavities with dominating Poynting flux (because weak fields would simply be anchored by the charges, 
hence would not decay) are familiar from the solar surface, and are correspondingly expected near the 
inner edges of stellar accretion disks, and near the innermost, strongly shearing galactic disks, in scaled- 
up proportions (Kundt, 1996, 2002). Such coronal magnetospheric reconnections in non-rigid rotators 
are expected to involve comparable powers to thermal emissions - cf. the magnetoid model of Ozernoy 
& Usov (1977) - because they drain on comparable energy reservoirs (controlled by equipartition). They 
can easily fill up the local hot bubble seen as the Broad Line Region (BLR) in AGN sources, which 
discharges to both sides of the disk through Blandford and Rees's (1974) deLaval nozzles, in the form of 
a supersonic twin-jet (Kundt, 1996), see Fig. 2. 

Once we deal with a central hot bubble filled with relativistic pair plasma, what will its energy dis- 
tribution be like? In the laboratory, atomic beams are routinely cooled by shooting at them with a laser 
beam. In the BLR, the charges are post- accelerated by the simultaneously generated Low-Frequency 
(LF) waves of the central magnetized rotator, of angular frequency uj reaching down to some 10~^s~\ 
hence of (large) strength parameter f : 

/ := eB/m^ c u = lO^^-^Bs/uj^i (10) 

for a typical coronal field strength B measured in KG; with B-s := B/IO'^G, u;_4 := u/lO^'^s^^. From 
the windzones of pulsars (like the Crab) we believe to have learned that the LF waves sweep the charges 
up in energy to Lorentz factors 7 of order p^^ = 10^-^(i?3/a;„4)^/'^, in the absence of damping (Kulsrud 
etal, 1972; Kundt, 1986). 

Damping occurs in the BLR through the equally present thermal (HF) radiation, a narrow 'bump' 
between IR and X-rays, via inverse-Compton losses which truncate a distribution towards high energies, 
because they scale as 7^. For class (d) (of AGN), these inverse-Compton losses are so strong that only 
some 10% of them are radio-loud, and show jets. 

As the outgoing charges in the (Thomson-opaque) BLR interact with both the LF and HF photons, 
their momentum distribution (away from the central engine) is expected to evolve towards a relativistic 
Maxwellian in radial direction, or even sharper, of Lorentz factor 7 > 10^, (corresponding to brightness 
temperatures in excess of lO^^K). Such high Lorentz factors are (i) expected, after equ. (10), are (ii) 
indicated by the overall energetics, (iii) by an avoidance of the inverse-Compton catastrophe, (iv) by the 
statistics of superluminal speeds (Kundt, 2004), and (v) by LF intraday variability, cf. Singal & Gopal- 
Krishna (1985), Wagner & Witzel (1995), but also Jauncey et al (2003) who prefer an interpretation via 
interstellar scintillations. 
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Figure 2: Sketch of a plausible Central Engine: Coronal Magnetic Reconnections create e''^ -pairs, the 
warm Central (Star and/or) Disk emits photons, Low-Frequency Waves post-accelerate the escaping e^, 
and the latter boost the thermal photons to high-energy 7-rays. An ambient thermal bulge serves as the 
deLaval nozzle from which a twin jet emerges, along the spin axis of the central rotator; in galactic-center 
sources, this region is observed as the BLR. 

When the charges leave the BLR and approach the deLaval nozzle - formed naturally by the obstruct- 
ing plasma of the ambient circumstellar medium, or central galactic bulge respectively, whose inertia (in 
pressure balance) scales inversely as its temperature T, i.e. is 10'^'^/T4-times larger - their narrow mo- 
mentum distribution will be channeled into a monoenergetic one, as inferred above from equ.(3), with 
dominating Poynting flux which carries at least half the energy. 



4. Discharging the Jets 

Once the highly-relativistic e^*" -pairs from the BLR enter a vacuum channel, rammed by preceding 
generations of charges, they form an almost loss-free, monoenergetic E x B-drifting beam as calculated 
explicitly in section 2 whose only losses are inverse-Compton collisions on the radiation background 
obeying: 



Ideg ■■= 7/7' = SnieC'^/iaTUsKl = MpC / 76(1 + z)^ (11) 

where Ideg is the degradation e-folding length, ar the Thomson cross section, u^k the energy density 
of the 3K background radiation, and 7e := 7/10^ , (Kundt, 2001). For a bulk Lorentz factor 7 ^ iC, 
thought to be realistic for most jet sources, inverse-Compton losses on Mpc scales are therefore ignorably 
small at redshifts z < 2. 

Note that a monoenergetic beam has no coUisional losses between its member charges - because 
they have vanishing relative velocities - nor dynamic-friction losses, for the same reason. Such internal- 
friction losses would in any case be ignorably small, because of the beam's extremely small electron- 
number density: 

Ue = L/A-fnieC^ = 10-Scm-3(L/A7)„3.5 . (12) 

During its propagation through a vacuum channel, a monoenergetic beam has only the inverse-Compton 
losses described by equ. (11). 

Conditions change when a beam encounters obstacles, in the form of (heavy) channel-wall material, 
or channel intruders, or obstructing material at its downstream end, its 'head'. Such obstructing plasma 
tends to be highly conductive, hence forbids penetration of electric and magnetic fields. The guiding 
toroidal magnetic field then gets compressed like the windings of a coil, and so are the convected charge 
clouds. Both electric charges and currents pile up against such a conducting wall, changing the field 
geometry in a way conveniently to be described by mirror charges (of opposite sign), and mirror magnetic 
fields (of same sign), see Fig. 3. 

Note that a charge-symmetric beam cannot be arrested, or reflected by electric fields alone; the latter 
can only redistribute the energies among the charges. The reflection of a (neutral) beam at its head is 
achieved by a changing geometry of both fields, electric and magnetic, such that the formerly quasi- 
stationary E X B-drift in forward direction is diverted, partially towards the beam's axis, and partially 
sideways towards its periphery (Fig. 3). At the same time, the bulk speed of the charges is reduced, 
from supersonic to subsonic, whereby straight-line motions change into gyrations. But now the charges 
are post- accelerated by the huge, convected electric potential, in the form of a space-charge limited flow 
whose relativistic version was first treated by Michel (1974), in application to pulsar polar-cap discharges. 

Michel's derivation of the relativistic generalization of Child's Law restricts itself to stationary, one- 
dimensional, one-fluid discharges inside a low-density plasma whose asymptotic speed is extremely rel- 
ativistic, and whose asymptotic charge density realizes force-freeness, also known as Goldreich- Julian 
(1969) density, or Hones-Bergeson (1965) density. The derivation takes care of the fact that due to dis- 
tributed screening, individual charges fall only through a tiny fraction of the available potential $. This 
fraction, of order the square root of e$ (in units of the electron rest energy), results as a consequence of 
Maxwell's equations plus conservation of energy, by integration along particle orbits through the (mag- 
netized) polar gap: 

7oo ~ (87$)'/' . (13) 
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Figure 3: Simplified cross section through a beam head, sketching the distributions of relativistic elec- 
trons (e^), electric and magnetic fields, mirror charges and fields (on the side of the obstructing ambient 
plasma), and particle orbits. Omitted are the motions of the stalled charges escaping from the impact 
center, whereby they are post-accelerated by the huge convected potential. 



The electron Lorentz factor 700 is reached exponentially towards the boundary of the polar gap, whose 
voltage $ is assumed stabilized by unipolar induction. 

In the present case of a stalled beam, the convected potential $ has its peak near the impact center, 
somewhat enhanced by compression w.r.t. its convected value, together with a surrounding ring wall 
of opposite sign. The boosting arena is therefore neither 1-d, 1 -fluid, nor stationary. Still, distributed 
screening will result in the square of the charges' Lorentz factor 7 tending (as above) towards some 
multiple of the maximum available one, 7$ , and we expect the energy distribution in the stalled beam to 
acquire a high-energy power-law tail - after averaging over a spatial ensemble of discharges - reaching 
up from its convected value (^beam ^ 10^) all the way to its peak value, of order (7$)^^^ ^ 10'' , see 
equ. (8). At the same time, the charges of inappropriate sign (and same instantaneous flow direction) 
are decelerated to lower energies, extending the spectrum downward in energy below 'jbeam- This sudden 



change of the energy distribution, from delta-like to hard power-law, takes place wherever a beam is 
stalled by ambient plasma, thanks to its convected Poynting flux. 

In this way, an almost loss-free, monoenergetic beam gets radiative whenever obstructed, with a broad 
power-law spectrum. A broad power law need not form, however, for peripheral tapping, in boundary- 
layer interactions. 

A well-known phenomenological dichotomy among the jet sources relates to the location of their 
hotspots, and tends to be called by their Fanaroff and Riley class I or II. Eilek et al (2002) have recently 
revised this classification into type A and type B, whereby class II is a subcase of type A, and speak of 
"straight" and "tailed" sources in the two cases. All jet sources start out "straight", their heads ramming 
supersonically into their ambient medium. During growth, the head's ram pressure drops as r~^ with 
distance r from the central engine, and the sound speed of the ambient medium often rises, so that beyond 
a certain distance - which differs from source to source - the head speed must pass from supersonic to 
subsonic w.r.t. the medium. From then on, the charges entering the terminating hotspot are no longer 
reflected (into their lobe) but continue coasting subsonically, in a gyrating mode, and can form a long, 
radiative "tail" (Gopal-Krishna et al, 1988, 1996). 

In our understanding, this sonic transition marks the transition from Eilek type A to B. Beyond the 
decelerating hotspot, particle motions are no longer channelled, or lossfree. The stalled jet material 
"engulfs" the ambient medium during relaxation, and squeezes it into small-filling-factor filaments, of 
huge relative mass densities pj : 

pJpH = QkT/niH^ = lO-'-^T-r (14) 

(in pressure balance, where X-ray temperatures T have been inserted, via T7). In this process, the pair 
plasma loses 1/3 of its injected energy, i.e. decelerates significantly, and "entrains" the thermal inclusions 
at a maximum speed of cpe/ pu = 10~^cTy , i.e. leaves them practically in their former state of motion. 
There is no beam beyond the terminating hotspot, yet there is ordered streaming at less than (2/3)c. 

5. Realistic Jets 

Our above treatment of astrophysical jets - being scale-invariant, and involving very similar central 
engines - is meant to apply to all (hundreds of) known jet sources: from (a) newly forming stars, (b) 
forming white dwarfs, (c) young neutron stars and BHCs (which latter are held to be neutron stars inside 
of heavy accretion disks), and from (d) the compact centers of (active) galactic disks; as elaborated in 
Kundt (1996, 2002, 2004). 

In all these cases, the central engine is (thought to be) a rotating magnet involving strongly sheared 
toroidal magnetic fields, whose intermittent reconnections supply abundant relativistic e'''-pair plasma, 
whose simultaneously emitted LFWs post-accelerate the pairs, and whose equally present thermal radi- 
ation cools the pairs, towards an outgoing relativistic Maxwellian and beyond. A deLaval nozzle forms 
automatically via the huge inertia of the ambient plasma, and the comoving Poynting flux sharpens the 



particle distribution to a monoenergetic one, and guides the pairs into the channels rammed by earlier 
generations. Such extremely relativistic beams are stabilized by their comoving Poynting flux, and by 
the inertia of the ambient plasma. This (generalized) 'unified scheme' for all jet sources deviates in part 
from most other approaches but compares favourably by being more explicit, more uniform, and more 
stable, and by not invoking (energy-rising) stochastic accelerations. In the sequel, we will comment on a 
few well- studied sources standing out by their extreme properties. 

The nearest active galactic nucleus is the center of our own galaxy, Sgr A*, the rotation center of 
the Milky Way disk. Its (unresolved) mass has been determined from the orbits of a few innermost 
stars as lO^'^^^'^Af©, its spectrum ranges from lO'^-^Hz to beyond 10 TeV, with an integrated luminosity 
of IO^^Lq peaking near GeV energies, and its output shows flares at IR and X-rays with variability 
timescales of ^ 17 min (Melia & Falcke, 2001; Genzel et al, 2003; Mayer-Has sel wander et al, 1998; 
Aharonian et al, 2004; Kundt, 1990, 2004; Roy & Rao, 2004). Are we seeing the innermost Galactic 
disk almost edge-on, whose present-day output reaches us at only 10^^ '^-times its Eddington value? Is its 
feeding presently throttled by fountain-like evaporation of the innermost disk? In any case, the radio part 
of the spectrum of Sgr A*, at z/ ^ lO^^Hz, of slope a := dlogSi,/dlogu = 0.3, signals monoenergetic 
synchrotron emission at Lorentz factor 7 ^ 10'', as does the Arc region (Anantharamaiah et al, 1991), 
whereas its enclosing emitter - Sgr A East - radiates a hard power-law spectrum at radio frequencies, as 
expected for stalled populations. 

An extreme case among extragalactic jet sources is the giant radio galaxy 3C 445, as concerns lossfree 
transport through large distances. Prieto et al (2002) emphasized the need for in-situ electron acceleration 
inside the hotspots, based on optical (synchrotron) emission nearly 0.3 Mpc away from the nucleus. From 
our equ. (11) it is clear that such a problem does not arise for E x B-drifting beams. See also Brunetti 
(2002), Hardcastle et al (2003), and Stawarz (2004) for similar well-studied sources. 

An even more extreme jet source is the quasar 3C 273, almost unique among 10'^ known radio jets by 
its brightness and one-sidedness (Morrison et al, 1984; Kundt & Gopal-Krishna, 1986; Jester et al, 2001, 
2002). Does its head plough almost luminally into its CGM (Ph ^ 0.6), and nearly towards us, with 
very little resistance in the (cosmic-ray?) halo of its host galaxy? Does its emitted spectrum soften on 
approach of its tip because of accumulating radiation losses, or because of increased beaming in forward 
direction (so that we observe a non-representative spectrum)? A rare source may well require a rare 
explanation. 

Proceeding to stellar jet sources, there is the unique binary neutron star, or BHC SS 433, whose 
interpretation has been controversial ever since its discovery as a jet source, in 1978 (cf. Fender, 2003; 
Kundt, 1996). Does it emit 'bullets' of local-galactic composition, at largely super-Eddington power, 
or are its beams composed of pair plasma, like in all the other jet sources? Kundt (2004) interprets its 
moving optical and X-ray emission lines as emitted by the impacted wind matter of its massive companion 
star, which forms its channel walls, and which is dragged along by the relativistic flow at a fraction of 
the speed of light. In the latter interpretation, we deal with a iC'-year young binary Galactic neutron 
star inside its SNR W 50, at a distance of 3 Kpc, whose (sub-Eddington) spindown power is still strong 
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enough to prevent accretion from its disk onto its surface, and to blow pairplasma jets whose heads have 
already crossed the periphery of W 50. 

That jets from Galactic binary neutron stars, or BHCs may consist of pairplasma has been recently 
advocated by Kaiser & Hannikainen (2002), via detections of the redshifted 511 KeV pair-annihilation 
line in eight binary X-ray sources with jets. When compared with that same emission line from the Crab 
pulsar (Massaro et al, 1991), its redshift (of some 7%) is more likely gravitational redshift from a neutron 
star's surface, where the density of slow pairs should be vastly higher than anywhere downstream along 
the jet. Another indication of highly relativistic electrons in compact stellar sources is the superluminal 
X-ray jet in the microquasar XTE J 1550-564, whose observed 'deceleration' may have to be understood 
as a varying phase velocity (Corbel et al, 2002). 

Finally, there are classes (b) (of PNe: Kundt, 1996; Balick & Frank, 2002), and (a) (of YSOs) which 
are hard to analyse because enshrouded by dense ionized, atomic, and/or molecular gas and dust. Among 
the few convincing (nonthermal) YSO candidates are two triple radio sources of (expansion) age ^ lO'^yr, 
one of them S 68 in Serpens (Rodriguez et al, 1989), further some 23 often one-sided core radio jets in 
stellar bipolar flows, including L 1455 (Schwartz et al, 1985), and HH 1 1 1 (Reipurth & Heathcote, 1993; 
Rodriguez & Reipurth, 1994; Reipurth & Bally, 2001), and the two-sided synchrotron jet from W3(0H) 
(Wilner et al, 1999). Note that in view of the opacity effects, a reliable distinction between thermal and 
synchrotron radiation has not always been possible. 
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